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1 Introduction 



A major phenomenon that the theory of strong interactions, quantum chromodynamics 
(QCD), should account for is confinement: quarks and gluons are not observable parti- 
cles. In fact, every physical process involving strong interactions at high energy results 
in the formation of hadrons, in which quarks and gluons are confined on a distance scale 
of 0(1) fm. While, up to now, there is no formal proof that QCD implies confinement, 
there are many indications, both from perturbative and from lattice numerical studies, 
that this is likely the case. Perturbative QCD is applicable to scattering processes of 
quarks and gluons involving large momentum transfer (^> 1 GeV) because the strong 
coupling constant as is small enough to allow a series expansion. However, this is no 
longer possible at a scale of 1 GeV or below, where the perturbative expansion is mean- 
ingless, and where confinement and hadronization, the process of hadron formation, 
takes place. Thus, hadronization is not yet calculable from QCD first principles and 
one has to resort to phenomenological models. While this may seem an inconvenient 
limitation, still much can be learned from these models about QCD in the confinement 
regime. Indeed, if they are able to effectively describe the essential features of the ac- 
tual physical process, they give us relevant information about the characteristics of the 
fundamental theory. 

In these lectures, we will review a model with a rather long history, that has re- 
cently been revived by its successes in the description of hadronic multiplicities both 
in high energy elementary collisions and relativistic heavy ion collisions: the Statistical 
Hadronization Model (SHM). 

2 The statistical hadronization model 

The idea of applying statistical concepts to the problem of multi-particle production 
in high energy collisions dates back to a work of Fermi [1] in 1950, who assumed that 
particles originated from an excited region evenly occupying all available phase space 
states. This was one of Fermi's favorite ideas and soon led to an intense effort in trying 
to work out the predictions of inclusive particle rates calculating, analytically and nu- 
merically, the involved multidimensional phase-space integrals. When it became clear 
that the (quasi) isotropic particle emission in the center-of-mass frame predicted by 
Fermi's model was ruled out by the data, an amendment was put forward by Hage- 
dorn [2] in the '60s, who postulated the existence of two hadron emitting sources flying 
apart longitudinally in the center-of-mass frame of a pp collision. Thereby, one could 
explain the striking difference between spectra in transverse and longitudinal momen- 
tum. Hagedorn was also able to explain the almost universal slope of pt spectra in his 
renowned statistical bootstrap model, assuming that resonances are made of hadrons 
and resonances in turn. 

After QCD turned up, many phenomenological models of strong interactions were no 
longer pursued and the statistical model was no exception. The resurgence of interest 
in these ideas came about when it was argued that a completely equilibrated hadron 
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gas would be a clear signature of the formation of a transient Quark-Gluon Plasma 
(QGP) in heavy ion collisions at high energy. While it has been indeed confirmed that 
an (almost) fully equilibrated hadron gas has been produced [3] in those collisions, the 
interest in this model was also revived by the unexpected observation that it is able 
to accurately reproduce particle multiplicities in elementary collisions [I]. Naively, one 
did not expect a statistical approach to work in an environment where the number of 
particles is O(10) because it was a belief of many that a hadronic thermalization process 
would take a long time if driven by hadronic collisions. Apparently this is not the case 
and one of the burning questions, which is still waiting a generally accepted answer, is 
why a supposedly non-thermal system exhibits a striking thermal behavior. 



Figure 1: High energy collisions are assumed to give rise to multiple clusters at the 
hadronization stage [top]. Each cluster [bottom] is a colorless extended massive ob- 
ject endowed with abelian charges (electric, strange, baryonic etc.), intrinsic angular 
momentum and other quantum numbers such as parity, C-parity and isospin. 

Before we address this interesting issue, it is appropriate to provide a rigorous for- 
mulation of the model in a modern form, which is necessarily different from Fermi's 
original model due to the tremendous improvement in our knowledge of strong interac- 
tions phenomenology. The Statistical Hadronization Model (SHM) must be considered 
as an effective model describing the process of hadron formation in high energy collisions 
at energy (or distance) scales where perturbative QCD is no longer applicable. A high 
energy collision is thought of as a complex dynamical process, governed by QCD, which 





Q = (Q,B,S,...) 
M = mass 

J, X = spin, 3 rd comp. 

V = volume 

Isospin, parity. C-oaritv 
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eventually gives rise to the formation of extended massive colorless objects defined as 
clusters or fireballs (see Fig. []]). While the multiplicity, masses, momenta and charges 
of these objects are determined by this complex dynamical process, the SHM postulates 
that hadrons are formed from the decay of each cluster in a purely statistical fashion, 
that is: 

Every multihadronic state localized within the cluster and compatible with 
conservation laws is equally likely. 

This is the urprinzip of the SHM. The assumption of the eventual formation of 
massive colorless clusters is common for many hadronization models (e.g. the cluster 
model implemented in the Monte-Carlo code HERWIG [5]) based on the property of 
color preconfmement [Hj exhibited by perturbative QCD. The distinctive feature of the 
SHM is that clusters have a finite spacial size. This aspect of clusters as a relativistic 
massive extended objects coincides with that of a bag in the MIT bag model [7]. Indeed, 
the SHM can be considered as an effective model to calculate bag decays. 

The requirement of finite spacial extension is crucial. If the SHM is to be an effective 
representation of the QCD-driven dynamical hadronization process, this characteristic 
must be ultimately related to the QCD fundamental scale Aqcd- As we will see, the 
universal soft scale shows up in the approximately constant energy density at hadroniza- 
tion; in other words, the volume of clusters is in a constant ratio with their mass when 
hadronization takes place. It is also worth stressing here that there is clear, indepen- 
dent evidence of the finite size of hadronic sources in high energy collisions. Quantum 
interference effects in the production of identical particles, the so-called Bose-Einstein 
correlations or Hanbury Brown-Twiss second-order interference, is by now a firmly es- 
tablished phenomenon. This effect would simply be impossible without a finite volume. 

3 Localized states 

The basic postulate of the Statistical Hadronization Model asserts that every localized 
multihadronic state which is contained within a cluster and is compatible with conser- 
vation laws is equally likely. The word localized, implying a finite spacial size, plays 
a crucial role, as we have emphasized. Thus, before getting to the heart of the SHM 
formalism, it is necessary to pause and clarify the distinction between localized and 
asymptotic states. 

Such a difference is not an issue when the volume is sufficiently larger than the 
Compton wavelength of hadrons and it is disregarded in most applications where clusters 
supposedly meet this requirement (e.g. heavy ion collisions); yet, it is an important 
point at a fundamental level. Although in thermodynamics the focus is on the limit of 
infinite volumes, we must start from a finite volume and localized systems to introduce 
concepts like energy density, temperature etc. Furthermore, in the hadronic world, finite 
size effects must be diligently taken into account when the volume is comparable to the 
(third power of) the pion's Compton wavelength, ~ 1.4 fm. 
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Figure 2: The localized multi-hadronic states \hy) pertaining to the quantum field 
problem in a limited region. Asymptotic states |/) are the usual free states characterized 
by particle momenta and spin components. 



The difference between a localized and an asymptotic state is depicted in Fig. [2j 
For a single particle in a Non-Relativistic Quantum Mechanics (NRQM) framework, the 
conceptual difference is easier to grasp: a localized state is described by a wavefunction 
which vanishes outside the cluster's region whereas an asymptotic state is a wavefunction 
which is defined over the whole space (e.g. a plane or a spherical wave). In Quantum 
Field Theory (QFT), a localized state is a state of the Hilbert space defined by the 
localized problem, e.g. the problem of the quantum field in the cluster's finite region. 
For a free field, if we enforce fixed or periodic boundary conditions, such states are simply 
defined by integer occupation numbers for each allowed mode in the finite region, as is 
well known. For a multiparticle state of non-interacting hadrons, this state will be 
defined by all occupation numbers of the modes determined by the fields associated 
with the different species of hadrons, and we will simply denote it with \hy) (where h 
stands for "hadrons" and V stands for the finite region of volume V). On the other 
hand, an asymptotic state is a state of the Hilbert space defined by the quantum field 
operators over the whole space; for a free field, these are the familiar multi-particle free 
states defined by, e.g., momentum and polarization: 

\f) = \Pl,<7l,-'-,PN,0N) (1) 

There is a noteworthy and deep difference between non- relativist ic quantum mechanical 
and quantum field theoretical case. In the latter, particle number is not fixed and 
a localized state with multiplicity (defined as the sum of all occupation numbers) N 
does not necessarily correspond to an asymptotic state with particle multiplicity TVq 
Unlike in NRQM, a localized state with multiplicity N has non-vanishing projections 
over asymptotic states with different particle multiplicities. In symbols: 



v — «o,7v|0) + a ltN \l) + . . . + o>n^n\N) + 



(2) 



1 We note that the multiplicity of an asymptotic free state is the properly defined number of particles 
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with the obvious condition that a^N —>■ when V — > oo for i ^ N . In particular, the 
vacuum of the finite-region problem |0)y is different from the vacuum of the full-space 
problem |0), which is commonly known as the Casimir effect. With a straightforward 
mapping of the Hilbert space of the localized quantum field onto the full Hilbert space, 
it is possible to express destruction operators of the localized field as linear combinations 
of destruction and creation operators of the field defined over the whole space [8] . These 
relations would be sufficient to calculate the coefficients of the above equation, but this 
is not really needed in the SHM, as it will be soon clear. 



4 The formalism: basics 

Let us consider a cluster and assume first that it can be described as a mixture of states. 
Then, the basic postulate implies that the corresponding density matrix is a sum over 
all localized states projected onto the initial cluster's quantum numbers: 

poc^PilMWPi^PiPvPi (3) 

where \hy) are multi-hadronic localized states and P, is the projector onto the cluster's 
initial conserved quantities: energy-momentum, intrinsic angular momentum and its 
third component, parity and the generators of inner symmetries of strong interactions 

a 

The operator Pj can be formally defined as the projector onto an irreducible vector of 
the full symmetry group and worked out in a group theory framework [HI [TDJ . It can be 
factorized into a "kinematic" projector, associated to general space-time symmetry, and 
a projector for inner symmetries. For the space-time symmetry, the relevant group is the 
extended orthochronous Poincare group 10(1,3)^ and an irreducible state is defined by 
a four-momentum P, a spin J and its third component A and a discrete parity quantum 
number ir = ±1. Therefore: 

Pj Pp, J,\,ir ^inner (4) 

If the projector Pp,j,A is worked out in the cluster's rest frame where P = (M, 0), it 
further factorizes into the product of simpler projectors [SHU], i-e.: 

Pp,j^ = 5\P-P)PjJ-±f^ (5) 

where P is the four-momentum operator, Pj\ is a projector onto SU(2) irreducible states 
| J, A) and II is the space reflection operator. 

As clusters are color singlets by definition, the projector Pi nner involves flavor and 
baryon number conservation. In principle, the largest symmetry group one should con- 
sider is SU(3), plus three other U(l) groups for baryon number, charm and beauty 

2 Opcrators in the Hilbert space will be denoted with a hat. Exceptions to this rule are projectors, 
which will be written in serif font. i.e. P. 
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conservation. However, SU(3) symmetry is badly broken by the mass difference be- 
tween strange and up, down quarks, so it is customary to take a reduced SU(2)®U(1) 
where SU(2) is associated with isospin and U(l) with strangeness. The isospin SU(2) 
symmetry is explicitly broken as well, but the breaking term is small and can generally 
be neglected. However, most calculations in the past have replaced isospin SU(2) with 
another U(l) group for electric charge, so that the symmetry scheme, from an original 

SU(2)j SOS pj n ®U(l) s ^ ran (j ieness ®U(l)h ar y 0n reduces to ^ iX) char ge®^ \X) strangeness®^ iX)baryon- 

Altogether, Pi nner can be written as 

Pinner = Pi,I 3 PqPx W 

where / and 7 3 are isospin and its third component, Q = (Qi, . . . , Qm) is a vector of 
M integer abelian charges (baryon number, strangeness, etc.) and P x is the projector 
onto C-parity, which makes sense only if the system is completely neutral, i.e. 1 = 
and Q = 0; in this case, P x commutes with all other projectors. 

From the density matrix ([3]) the probability of observing an asymptotic multiparticle 
state \f) is 

P/OctflPiPyP,!/) (7) 

which is well-defined in terms of positivity and conservation laws. In fact, (J7j) is mani- 
festly positive definite and pf = if the state |/) has not the same quantum numbers as 
the initial state. By summing over all states |/), one obtains the trace of the operator 
PiPyPj which is 

J^Pf °c tr(PiP v Pi) = tr(P, 2 PvO = atr(P*P y ) . (8) 

/ 

The constant a is divergent and positive. It can be directly checked by choosing the |/) 
as momentum eigenstates and using the expression on the right hand side of (j3j). The 
reason for its presence is the non-compactness of the Poincare group, which makes it 
impossible to have a properly normalized projector. The last trace in (jSJ) can be written 

as 

tT(p i p v ) = ^2{h v \p i \h v ) = n (9) 

h v 

which is, by definition the microcanonical partition function [5], i.e. the sum over all 
localized states projected onto the conserved quantities defined by the selected initial 
state. If only energy and momentum conservation is enforced, Q takes on a more familiar 
form: 

tt = ^2(h v \5 4 (P-P)\h v ). (10) 

hv 

Although the mixture of states defined by Eq. ([3]) allows us to calculate probabilities 
of any measurement unambiguously, a cluster could in principle also be described by a 
pure quantum state. Actually, the mixture of states only expresses our ignorance about 
the true state of the system, which is in principle a pure one, or, more precisely, a pure 
state entangled with pure states of other clusters. To avoid slipping into fundamental 
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quantum mechanics problems of decoherence and measurement, we take a pragmatic 
stance here. It suffices to realize that in some low-energy collision events, only one 
cluster might be created whose state is then necessarily a pure one. According to the 
postulate of the SHM, this must be an even superposition of all localized states with 
the initial conserved quantities, i.e. 

\ip) = Ch v Pj\hy) with \ch v \ 2 = const . (11) 

h v 

The probability of observing a final state |/) is then 

iw>i 2 = iE</i p «i^vi a (12) 

hv 

= const X)K/l p <IMI a + Yl (f\ P i\ h v)(hv\Pi\f)ch v c* hlv . 

hy hy^h'y 

If the coefficients Ch v have random phases, the last term in Eq. (JT2J) vanishes and we are 
left with the same expression appearing in Eq. ([7]); in other words an effective mixture 
description is recovered. Hence a new hypothesis is introduced in the SHM here: if the 
cluster is a pure state, the superposition of multi-hadronic localized states must have 
random phases. 

Now the main goal of the model is to determine the probabilities flTJ) which involves 
the calculation of the projector Py = J2h v \hy){hy\, a more limited task than the 
explicit calculation of all scalar products (hy\f). Since the states \hy) are a complete 
set of states of the Hilbert space Hy for the localized problem, the above projector is 
simply a resolution of the identity of the localized problem and can be written in the 
basis of the field states. For a real scalar field this is 

Py= [ (13) 

where \ip) = ® X |V ; ( X )) an d Dip is the functional measure; the index V means that the 
functional integration must be performed over the field degrees of freedom in the region 
V, that is Vip = Ylxev dip(~x). One has to face several conceptual subtleties in the 
endeavor of calculating the probabilities ((7j) with the projector (fT5|) . e.g. how to deal 
with field boundary conditions and with their values outside the region V. However, by 
enforcing the known non-relativistic limit is possible to come to an unambiguous and 
consistent result [8]. 

5 Rates of multiparticle channels 

According to the formulae introduced in the previous section, the decay rate of a cluster 
into a channel {Nj} ({Nj} is the array of multiplicities (Ni, . . . , Nk) for hadron species 
1,...,K) is proportional to the right hand side of (|7j) integrated over momenta and 
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summed over polarization states of the final hadrons. Taking into account only energy- 
momentum conservation, so that the projector (jSj) reduces to 

Ppjx, ^Pp = S 4 (P - P) , 

and neglecting quantum statistics effects, this is proportional to the microcanonical 
partition function with fixed particle multiplicities [91 [8] 



n 



v 



N 



{Nj} 



(2S j 



(2k) 



3iV 



J d 3 Pl ...J d 3 PN S\P -Y,P*)(U\Pv\0). (14) 



Here N is the number of particles, Sj the spin, Pq = (M, 0), M is the mass and V is the 
cluster's proper volume. This formula is the same as it would be obtained in NRQM, 
with the factor (0|Py|0) (which becomes 1 in the limit V — > oo) being the only effect 
of the field theoretical treatment [Sj. Since only relative rates make sense, this common 
factor for all channels is irrelevant. 

Loosely speaking, Eq. (fT4l) tells us that the decay rate of a massive cluster into some 
multi-hadronic channel is proportional to its phase space volume. However, it should 
be emphasized that the "phase space volume" in (1141) is calculated with the measure 
d 3 xd 3 p for each particle, and not with the one usually understood in QFT, i.e. d 3 p/2e. 
Although this is also commonly known as "phase space", it is quantitatively different 
from the properly called phase space measure d 3 x d 3 p and should be called " invariant 
momentum space" measure |12j . 

Eq. f[T4"|) can be cast in a form which makes its Lorentz invariance apparent. Define 
a four-volume T = Vu [12] where V is the cluster's rest frame and u its four-velocity 
vector. Then ffT4D can be rewritten as: 



n 




(2tt) 3JV 

N 



d A p x ... \ d A p N 



(15) 



5 4 (P -$>)(0|Py|0> 



which is manifestly covariant. In this form it can be directly compared with the general 
formula for the decay rate of a massive particle into a iV-body channel: 



Tat oc ^ 




d 3 p N 



2s 



\M 



N 



fi 



2 *\Po - E 



Pi 



(16) 



where a labels, as usual, polarization states. Comparing (fl4"|) with fflBI) we can infer a 
dynamical matrix element for the SHM which is 



N 



\M 



fi I 



X 



nx 

i=l 



■Pi 



(17) 
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Therefore, according to the SHM the dynamics in cluster decay is limited to a common 
factor for each emitted particle, which linearly depends on the cluster's spacial size. The 
four-volume T is simply proportional to the four-momentum of the cluster through the 
inverse of energy density p and therefore: 



1 N 

jvlp-**" ( 18 ) 



\M H \ 2 oc 
I fi\ p 



This expression explicitly shows the separation between the kinematic arguments of the 
dynamical matrix element, and the scale 1/p which determines particle production. This 
ought to be ultimately related to the fundamental scale of quantum chromodynamics, 
Aqcd- 

It has already been stressed that the finite cluster size is the distinctive feature of 
the SHM. This peculiarity of the model stands out when taking into account quantum 
statistics for the calculation of decay rates. Our final result, for which (TTJJ is a special 
case when all particles belong to different species, reads 

<w = / *pi . . . <*w s\ Po - n e {T1) Z +Hj{2 " j + (19) 

■' < i 3 {h nj } II,.;' i », ■ 

Hi 

[[F (0|Pv|0) 



X 



L=l 



where {h n } is a partition of the integer Nj in the multiplicity representation, that is 
Fourier integrals: 



Yln'=i n jh nj = Nj, Y^n-=ih nj = Hj and E 7 ^ = N - Tlie factor F n; . in Eq. (JT9D are 



1 r 

F «« = 11(2^3 jf ^e^W^ (20) 

over the cluster's region V, c\ being a cyclic permutation of order n\. The expression 
f fT9|) has been obtained in refs. [HI [8] and is a generalization of a similar one calculated 
by Chaichian, Hagedorn and Hayashi [12] whose validity is restricted to large volumes. 
It is a so-called cluster decomposition of the microcanonical partition function of the 
channel. For sufficiently large volumes, all terms in Eq. fTTTJT) turn out to be proportional 
to the Hjth power of the volume V [9] , so that the leading term is the one with Hj = Nj 
for all j, which leads precisely to Eq. (fT4l) . Thus, Eq. (fT9l) is a generalization of <HM 
containing all corrections due to quantum statistics. 

In general, with respect to the Boltzmann case (j!4p . the channel rate is enhanced 
by the presence of identical bosons and suppressed by that of fermions. This means 
that Bose-Einstein and Fermi-Dirac correlations are built into the SHM. The reader has 
probably anticipated this fact through the appearance of typical Fourier integrals in the 
cluster decomposition. This feature of the SHM is an almost obvious consequence of the 
cluster's finite spacial size. 
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6 Interactions 

So far, we have dealt with non-interacting particles. However, the localized hadronic 
fields we have used to calculate transition probabilities do interact and this must be 
taken into account. The energy of the interacting system must be conserved until the 
final asymptotic multi-hadronic state is reached which is made of particles stable under 
strong interactions, namely pions, kaons, nucleons and octet hyperons. 

Formally, this implies that the projector ([5]) must include the interacting Hamil- 
tonian in the <5 4 (P — P) operatoiiE The definition (EJ) for the probability to observe 
an asymptotic state |/), has to be modified by the insertion of M0ller operator fl, yet 
summing over the complete set of states yields the same result as in Eq. (jSJ): 



where a is an irrelevant divergent constant and f2 is the microcanonical partition function 
of the interacting hadronic system. 

An outstanding theorem by Dashen, Ma and Bernstein (DMB) [Hj allows us to 
calculate the microcanonical partition function of an interacting system in the thermo- 
dynamic limit V — > oo as the sum of the free one plus a term depending only on the 
physical scattering matrix. It can be expressed as 



where P includes the full interaction Hamiltonian, whereas Pq only contains the free one; 
S is the reduced scattering matrix on the energy-momentum shell. If more conserved 
quantities other than energy and momentum are involved, like those encountered in 
section [H the theorem is readily extended and relevant projectors can be placed next 
to the ^-functions in Eq. (1221) ; it suffices that these conserved quantities are associated 
with symmetries of both free and interacting theory. 

This theorem is indeed the starting point of the hadron-resonance gas model since 
it can be shown that if only the resonant part of the scattering matrix is retained and 
the background interaction can be neglected, the main contribution of the second term 
on the right hand side of Eq. (1221) is equivalent to considering all hadronic resonances 
as free particles with distributed mass. More specifically, if a cluster decomposition of 
the two scattering operators is carried out in Eq. (j22|) . the corresponding diagrams can 
be divided in two sets: the symmetric diagrams (see Fig. [31 left panel) and the non- 
symmetric ones (see Fig. [31 right panel). Taking into account that the terminal legs on 
both sides have to be the same stable particles on entry and exit (we are calculating 
a trace), it can be shown that the main contribution to symmetric diagrams comes 

3 In all virtually known field theories, there is no additional interacting term for momentum and 
angular momentum |13| . 






(22) 
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Figure 3: Left panel: symmetric diagrams for the cluster decomposition of the interaction 
term in the DMB theorem. Right panel: non-symmetric diagrams. 

from the matching resonances in bubbles facing each other. For each term of the trace, 
this amounts to adding the decay products of resonances considered as free particles 
with masses distributed according to a relativistic Breit-Wigner form. In symmetric 
diagrams, there is in principle an additional contribution from resonance interference, 
which might be non-negligible in case of wide, overlapping resonances with the same 
decay channel, but it depends on mostly unknown complex parameters and it is thus 
neglected. 

Likewise, the asymmetric diagrams give an additional contribution which also de- 
pends on the aforementioned complex interference parameters. While the number of 
such diagrams greatly exceeds the symmetric ones due to the large number of reso- 
nances, contributing terms can be both positive and negative and hopefully a partial 
cancellation occurs when summing them up for a selected final state. 

Altogether, retaining only the resonant interaction and symmetric diagrams in the 
scattering matrix cluster decomposition, and neglecting resonance interference leads to 
the following picture: an interacting hadron gas is, to a good approximation, a gas of 
non-interacting free hadrons and resonances. Since non-resonant interaction should be 
negligible, the ideal hadron-resonance gas picture holds if the energy density or temper- 
ature of the system is large enough for most resonances to be excited. A quantitative 
assessment of how large these parameters are is still missing, a rough estimate being 
T > 100 MeV. 

An important remark is now in order. The DMB theorem affirms the equality of two 
traces, but not of single trace terms. Yet, the decomposition of Eq. (1221) . implying the 
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ideal hadron-resonance gas picture, is widely used for the calculation of inclusive stable 
hadronic multiplicities as well, which requires a condition stronger than the equality 
of the traces on both sides. In other words, using the decomposition (122]) to calculate 
average multiplicities or fluctuations requires the equality to hold for multiparticle gen- 
erating functions and not only for microcanonical partition functions. Up to now, the 
extension of (1221) to generating functions has never been proved; most likely, it does not 
hold and corrections to this assumption are necessary. Moreover, while the theorem re- 
quires the thermodynamic limit, it is commonly used at finite volume. These limitations 
should be always kept in mind when using the ideal hadron-resonance gas model. 

7 High energy collisions 

As we have seen in Sect. O each individual cluster produced in a high energy collision 
(shown in Fig. [1]), should be hadronized according to formula ( fT9|) . or its approximation 
(TH|) . which yields the rates of a given channel within the microcanonical ensemble, in- 
cluding energy-momentum conservation. If clusters are large enough, the microcanonical 
ensemble could be well approximated by the canonical [9j [15] or even grand-canonical 
ensemble for average multiplicities. This is not the case in elementary collisions (e + e~, 
pp, etc.) while it is generally possible in heavy ion collisions, as we will see later in this 
section. Calculating observables in high energy collisions within the SHM then implies 
summing microcanonical averages over all produced clusters and this requires in turn 
knowledge of their charges and four-momenta. In fact, this latter information is un- 
known to the SHM and only a dynamical model of the pre-hadronization stage of the 
process (such as, e.g., HERWIG) can provide it. 

However, if we are interested in calculating Lorentz-invariant observables (such as 
average multiplicities) the momenta of clusters are immaterial and only charges and 
masses matter. In this case, one can introduce a peculiar extra-assumption which allows 
to considerably simplify the calculation. Basically, it is assumed that the probability 
distribution: 

w{Q 1 ,M 1 ;...,Q N ,M N ) 

of masses M and conserved abelian charges Q for N different clusters is the same 
as one would have by randomly splitting a large cluster (defined as Equivalent Global 
Cluster, EGC) into N subsystems with given volumes. Thereby, the Lorentz invariant 
observables can be calculated for one (equivalent global) cluster, whose volume is the 
sum of proper cluster volumes and whose charge is the sum of cluster charges, hence 
the conserved charge of the initial colliding system. The full mathematical procedure is 
described in detail in ref. [To] . 

In such a global averaging process, the EGC generally turns out to be large enough 
in mass and volume so that the canonical ensemble becomes a good approximation 
of the more fundamental microcanonical ensemble [15]; in other words, a temperature 
can be introduced which replaces the a priori more fundamental description in terms 
of energy density. It was shown that the mass of the EGC should be at least 8 GeV 
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(with an energy density of 0.5 GeV/fm 3 ) for the canonical ensemble to be a reasonably 
good approximation [15]. Also, it should be emphasized that in such a mathematical 
reduction process, temperature has essentially a global meaning and not local as in 
hydrodynamical models (see next subsection). The only meaningful local quantity in 
actual physical process are energy densities and individual physical clusters cannot be 
described in terms of a temperature, unless they are sufficiently large. Nevertheless, 
this "global" temperature closely mirrors the value of energy density at which clusters 
hadronize. Indeed, it is this latter value which mainly determines hadronization-related 
observables; the requirement of the charge distribution of EGC is a side-assumption 
which is important to simplify calculations, but it can possibly be replaced by other 
distributions leaving final results essentially unchanged. 

In this approach, the primary multiplicity of each hadron species j is given by 



= TO + i) | ;7? ,„ (Ti) „ +1 ^ Kj ^ (23) 



n=l 



where V is the (mean) volume and T the temperature of the equivalent global clus- 
ter. Here Z(Q) is the canonical partition function depending on the initial abelian 
charges Q = (Q, N, S,C, B), i.e., electric charge, baryon number, strangeness, charm 
and beauty, respectively; rrij and Sj are the mass and the spin of the hadron j, cy = 
(Qj, Nj, Sj, Cj, Bj) its corresponding charges; the upper sign applies to bosons and the 
lower sign to fermions. 

The parameter 75 in (|23|) is an extra phenomenological factor implementing an ad hoc 
suppression of hadrons with N s strange valence quarks with respect to the equilibrium 
value. This parameter is outside a pure thermodynamical framework and it is needed 
to reproduce the data, as we will see. For temperature values of 160 MeV or higher, 
Boltzmann statistics, corresponding to the term n = 1 only in the series (1231) . is a very 
good approximation (within 1.5%) for all hadrons but pions. For resonances, the formula 
( 123]) is folded with a relativistic Breit-Wigner distribution for the mass rrij. The final 
multiplicities, to be compared with the data, are determined by adding to the primary 
multiplicity (123]) the contribution from the decay of unstable heavier hadrons, according 
to the formula 

fa) = (n 4 )P~y + ^Br(A;^j)K) • (24) 
The canonical partition function can be expressed as a multi-dimensional integral 



x exp 



V 



(2tt)3 



^(25^ + 1) J d 3 p log (1 ± ^e-V^P^-^ 



±1 



(25) 



where N is the number of conserved abelian charges. Unlike the grand-canonical case, 
the logarithm of the canonical partition function does not scale linearly with the volume. 
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Figure 4: Upper panel: measured vs theoretical multiplicities of light-flavoured hadrons 
in e + e~collisions at ^ = 91.25 GeV. Lower panel: fit residuals (from ref. [T7]). 



Therefore, the so-called chemical factors Z (Q — nqj) / Z (Q) [18] turn out to be less than 
unity even for a completely neutral system at finite volume (canonical suppression) and 
reach their grand-canonical value 1 at asymptotically large volumes [HI H] . 

The light-flavoured multiplicities in e + e~show a very good agreement with the pre- 
dictions of the model, as it shown in Fig. HI the temperature value is about 160 MeV 
and the strangeness undersaturation parameter 75 ~ 0.7. Similar good agreements are 
found for many kinds of high energy elementary collisions over a large energy range [4]. 
Also, an excellent agreement between measured and predicted relative abundances of 
heavy flavoured hadronic species in e + e~collisions by using the model parameters fitted 
to light-flavoured multiplicities [HE], as shown in Tabled] 

The overall striking feature is that the temperature turns out to be approximately 
constant over two orders of magnitude in centre-of-mass energy with a value of 160- 
170 MeV (see Fig. [5]) and very close to the QCD critical temperature as determined 
from lattice calculations. There must certainly be a profound connection between the 
thus-found hadronization temperature and QCD thermodynamics, a connection which 
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Figure 5: Temperatures fitted in elementary collisions as a function of center-of-mass 
energy. 



has not been made clear yet. Nevertheless, this finding indicates that hadronization is 
a universal process occurring at a critical value of the local energy density, i.e. when 
clusters have an energy density of ~ 0.5 GeV/fm 3 . 

The parameter 75 is found to be less than 1 in all examined elementary collisions, 
ranging from ~ 0.5 in hadronic collisions to ~ 0.7 in e + e~ collisions (see Fig. [TOl) . This 
extra parameter most likely reflects the different mass of the strange quark with respect 
to lighter u, d quarks. This is a second scale, besides Aq C d, which must play a role 
in hadronization in view of its value 0(100) MeV. Altogether, one can say that the 
SHM description of hadronization is in excellent agreement with QCD at least with 
regard to the the number of parameters. The two parameters T and 75 correspond to 
the two fundamental scales Aqcd and m s , the strange quark mass. While we lack a 
definite relation connecting them (see, however ref. [T7j), it is worth stressing that a 
phenomenological description of hadronization in terms of fewer parameters cannot be 
possible. 

Finally, the statistical model shows a very good capability of reproducing transverse 
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Table 1: Abundances of charmed hadrons in e + e~ — > cc annihilations and bottomed 
hadrons in e + e~ — > bb annihilations at y/s = 91.25 GeV, compared to the prediction of 
the statistical model (from ref. [T7]). 



momentum spectra in hadronic [16] as well as heavy ion collisions [20] (see fig. [6]). 
Particularly, the phenomenon of approximate my scaling observed in pp collisions is 
nicely accounted for by the model. However, the exact conservation at low energy 
and the increasing importance of jet emission at high energy restrict the validity of 
the statistical canonical formulae to a limited centre-of-mass energy range. Within this 
region, a clear consistency is found between the temperature parameter extracted from 
the spectra and that from average multiplicities. Altogether, this finding bears out one 
of the key predictions of the SHM, namely the existence of a definite relation between 
the dependence of particle production rates on mass and, for each particle species, their 
momentum spectra (in the cluster's rest frame) because they are both governed by one 
parameter, the energy density (or temperature) at hadronization. 
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Figure 6: Transverse momentum spectrum of ir + in pp collisions at y/s = 27 GeV 
(from ref. [16]). Full dots are data from experiment NA27; the black line is a fit with 
temperature T=161 MeV and average transverse four- velocity of hadronizing clusters 
~ 0.21. Coloured lines show the cumulative contribution of resonance decays, divided 
into classes according to their quantum numbers. 



8 Heavy ion collisions 

In heavy ion collisions, the system is much larger and two possibilities are usually en- 
visaged: either hadronizing clusters are simply much larger than those in elementary 
collisions; or clusters are hydrodynamical cells, i.e. they are small but in thermal con- 
tact with each other due to previous thermalization, which implies a strong correlation 
between their position and momentum and charge densities (see Fig. [7]). In both case 
the canonical or grand-canonical formalisms apply to individual clusters. For the former 
case, it is worth mentioning that the transition from a canonical to a grand-canonical 
description effectively occurs when the cluster volume is of the order of 100 fm 3 at an 
energy density of 0.5 GeV/fm 3 [18]. If the EGC reduction assumption still applies, 
chemical factors in Eq. fT23|) are replaced by fugacities, and in this case the phase-space 
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Figure 7: Spacial distributions of clusters in heavy ion collisions according to the hy- 
drodynamical picture. In this model, nearby clusters interact from an early stage on 
and their momenta and charges are strongly correlated with their positions, unlike in 
elementary collisions. 



integrated multiplicities read 

K-> primary = VTi l% +1) E 7^ n (Tir +1 ^ K 2 (^) exp[n„ • q,/T] . (26) 



n=l 



fi is a vector of chemical potentials pertaining to the conserved abelian charges, i.e. the 
electrical chemical potential /xq, the baryon chemical potential [Ib and the strangeness 
chemical potential fis- Usually, but not always, /is and /iq are determined by enforcing 
strangeness neutrality and by fixing the ratio Q/B to be the same as the initial Z/A 
ratio of the colliding nuclei. 

In the framework of the hydrodynamical model, formula (1261) applies to individual 
clusters identified with hydrodynamic cells and both temperature and chemical poten- 
tials depend on space-time; when integrating particle densities to get average multiplic- 
ities, one should take into account this dependence. It is important to stress that the 
hydrodynamical description is a salient feature of heavy ion collisions due to the early 
thermalization of the system in the partonic phase, a phenomenon which does not oc- 
cur in elementary collisions. It is this early thermalization which establishes the strong 
correlation between positions and velocities of clusters, supposedly absent in elementary 
collisions. 

Provided that rapidity distributions are wide enough, and that there is little variation 
of the thermodynamical parameters of clusters around midrapidity, the formula (T2"6"D 
describes rapidity densities of hadrons at midrapidity as well: this condition is fulfilled 
at RHIC energies, but not at AGS and SPS energies, where the measured rapidity 
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Figure 8: Upper panel: measured vs theoretical multiplicities of light-flavoured hadrons 
in Pb-Pb collisions at \^s NN = 17.2 GeV. Lower panel: fit residuals (from ref. |21j). 



distributions are not significantly wider than those of a single fireball at the temperature 
found [22]. 

In general, the fits to particle multiplicities in heavy ion collisions are of the same 
good quality as in elementary collisions (see Fig. [8]). Many groups have analyzed the 
data over more than a decade [3] and the overall description is very good throughout 
all explored energies and one finds a smooth curve in the T — fiB plane (see Fig. [9]) . 

9 Strangeness production 

The statistical model is a very useful tool to study one of the main features of relativis- 
tic heavy-ion collisions, the increase of relative strangeness production with respect to 
elementary collisions. This was one of the early signatures proposed for Quark-Gluon 
Plasma formation [21], and it has therefor attracted much attention both on the theoret- 
ical and experimental side. According to the SHM, this is mainly an effect of the increase 
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of the global volume from elementary to heavy-ion collisions. In elementary collisions, 
the EGC volume is small enough for the chemical factors (see Eq. ( 1251) ) of strange 
particles to be consistently less than 1 for systems with vanishing net strangeness, a 
phenomenon known as strangeness canonical suppression. 

However, canonical suppression is not enough to account for strangeness enhance- 
ment from pp to heavy-ion collisions: also an increase of 75 is needed. This is demon- 
strated by neutral mesons containing strange quarks, especially meson, which do not 
suffer canonical suppression but are relatively more abundant in heavy-ion collisions 
[251 126] . Therefore, from a SHM viewpoint, one can say that, as far as particle abun- 
dances is concerned, the only substantial difference between elementary and heavy ion 
collisions resides in the different 75 values, which are generally higher in heavy-ion colli- 
sions and increase slowly as a function of center-of-mass energy (see Fig. [TU]) : at RHIC 
energies one finds 75 ~ 1 in central collisions. However, since 75 is an empirical param- 
eter which lies outside of a pure statistical mechanics framework, this observation does 
not clarify the origin of strangeness enhancement. 

It is interesting to note that, while 75 shows no special regularity in elementary 
collisions, the ratio of newly produced ss pairs over one half uu + dd pairs (the so- 
called Wroblewski ratio As) turns out to be around 0.2-0.25 at all energies, whereas it 
is definitely higher in heavy ion collisions (see Fig. ITTT) . 

These differences, both in 75 and Xs have been and still are subject of investigation. 
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of center-of-mass energy in e + e~, pp and pp collisions. Right panel: the strangeness 
undersaturation parameter 75 as a function of center-of-mass energy in central heavy 
ion collisions (from ref. [23]). 



Relevant information comes from the centrality dependence of strangeness production, 
which provides an interpolation from single pp collisions at large values of nuclear impact 
parameter to head-on heavy-ion collisions at low values thereof. The enhancement has 
been measured by the experiments WA97 and NA57 at SPS energy [27] and STAR 
at RHIC [28] for hyperons and other strange particles and it has been found to be 
hierarchical in strangeness content (highest for Q~ , lowest for A). These observations 
led some authors [2H] to put forward a picture where 75 is an effective parametrization 
of a canonical suppression. For large enough baryon number and charge, it is possible to 
take a mixed canonical-grand-canonical approach where only strangeness conservation 
is enforced, while electric and baryon-chemical potential are introduced. The chemical 
factors Z(S—Sj)/Z(S) depend on the volume and saturate at large volumes, as expected. 
Therefore, if we want to account for 75 < 1 with this mechanism, there must be some 
small sub-regions within a large fireball where strangeness is exactly vanishing even for 
the most central collisions. Thereby, chemical factors are significantly less than 1 and a 
suppression with respect to the grand-canonical limit is implied. However, this model 
has two major problems: 

1 . Since measured enhancement steadily increases from peripheral to central collisions 
and hadronization temperature does not change [301 EH E3], the volume of the 
sub-regions with S = should also increase and a saturation is thus expected (see 
Fig. [T2l : yet, no saturation is observed, which is quite an oddity. 

2. As has been mentioned, canonical suppression has no effect on <fr; yet, the relative 
yield of this meson with two constituent strange quarks is also observed to increase 
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Figure 11: Wroblewski ratio Xs (see text for definition) as determined in elementary 
and heavy ion collisions from fits of multiplicities to the statistical hadronization model. 



from peripheral to central collisions [32j and with 75 = 1 and the observed constant 
temperature, this cannot occur. 

Recently, a geometrical explanation of these two features has been advocated [33] 
based on a superposition of emission from a hadron-resonance gas at full chemical equi- 
librium with 75 = 1, defined as the core, and from nucleon-nucleon collisions at the 
boundary of the overlapping region of the two colliding nuclei, defined as the corona, 
from which produced particles escape unscathed. Since in NN collisions strangeness is 
suppressed with respect to a fully equilibrated, grand-canonical hadron gas, if such NN 
collisions account for a significant fraction of total particle production, a global fit to 
particle multiplicities will find 75 < 1, as indeed observed in data. The idea of superpos- 
ing different sources is common to other models (a similar one is discussed in ref. [34"]). 
The peculiar feature of this specific model is to assume single NN collisions as secondary 
sources; only in this case does it seem possible to reproduce centrality dependence of 
the (j) meson. 
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Figure 12: Canonical enhancement (defined here as the ratio between the chemical 
factor Z(S — Sj)/Z(S) and its value at some fixed volume Vo) as a function of volume 
for hyperons (from ref. [29]). 

10 Thermalization: how is it achieved? 

After discussing the success of the SHM in reproducing particle multiplicities and the 
intriguing universality of its main parameter, the temperature, one is obviously led to 
the question how this can come about. A classical process of thermalization through 
binary collisions between formed hadrons, advocated in heavy ion collisions [35], is 
ruled out in elementary collisions because the expansion rate is fast and hadrons are 
not interacting for a long enough time for this to take place. But even in heavy ion 
collisions peculiar features of the data cannot be explained in a hadronic kinetic picture 
[36] without invoking the predominance of multi-body collisions; since, in this case, 
the hadronic mean free path is comparable or smaller than Compton wavelengths, the 
collisional picture breaks down naturally. 

There is evidence that thermalization occurs at a relatively early stage over a large 
region (i.e. clusters several femtometers wide) in heavy ion collisions, whereas it is a 
late phenomenon (i.e. very close to hadronization) occurring over small (of the order of 
1 fm) distances in elementary collisions. Yet, the agreement between model and data 
is surprisingly accurate in elementary collisions, even more accurate than in heavy ion 
collisions, the only difference being in the level of strangeness phase space saturation. 
Somehow, the hadrons must be born into equilibrium as Hagedorn first pointed out [37] 
and was reaffirmed by others [381 EH] • 
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The idea that this thermal-like behavior is of genuine quantum-mechanical origin 
and not related to semi-classical collision processes, is shared by many [39] including the 
author. A different point of view was presented in ref. [10] where it was argued that the 
thermal behavior could just be mimickedby a matrix element which is weakly dependent 
on the final kinematic variables in (116j) . a scenario called "phase space dominance". 
But even this scenario requires stringent conditions on the dependence of cluster decay 
rates on the channel multiplicity (essentially like A ) |10] such that the exponential 
dependence of production rates on mass is not spoiled. Hence phase-space dominance is 
not less trivial in any way. A possible path to distinguish between the two scenarios is 
provided by the analysis of exclusive rates at low energy, although it must be pointed out 
that the observed identical particle correlations already favors SHM which is endowed 
with a built-in spacial extension, unlike phase-space dominance. 

However, whether it is a proper thermal-statistical equilibrium in a finite volume 
or rather a phase space dominance effect, there must be a profound reason behind this 
phenomenon, which ultimately has to be related to the nature of QCD as a theory with 
strong coupling at large distances. Also, we believe that the intriguing universality of 
the temperature found in elementary collisions as well as heavy ion collisions and its 
resemblance of the QCD critical temperature is not accidental. 

If we assume that post-hadronization collisions are unable to restore equilibrium, 
how can a quantum evolution process ensure it? Several years ago it was pointed out 
that a closed quantum system whose classical counterpart is chaotic and ergodic can 
give rise to thermal distributions provided that the so-called Berry conjecture applies 
|41j . Berry's conjecture [12] essentially states that the high- lying eigenfunction am- 
plitudes ^(x) in configuration space appear to be random Gaussian numbers and, as 
a consequence, momentum space distribution is microcanonical [13]. This "quantum 
thermalization" mechanism has been invoked to explain the observed thermal-like dis- 
tributions in hadronic processes [2]. Of course, this argument requires that classical 
QCD is chaotic (as it has indeed been advocated [15]), that Berry's conjecture holds for 
quantum fields, and that it can be applied to a dynamical process like hadronization. 
All of these conditions are non-trivial, but pursuing these ideas further may give rise to 
interesting developments. 

Recently, another appealing idea to explain the universality of thermal features in 
multihadron production has been put forward It invokes an analogy between con- 
finement and black hole physics. It is conjectured that the phenomenon of confinement 
is equivalent to the formation of an event horizon for colored signals (quarks and glu- 
ons). Similarly to Hawking- Unruh radiation, the spectrum of hadrons, emitted as the 
result of a high energy collision, is thermal because no information can be conveyed 
from the causally disconnected region beyond event horizon. According to this so-called 
Hawking-Unruh scheme of hadronization, temperature is related to the string tension 
and is thereby universal. Another interesting consequence of this idea is that the extra 
strangeness suppression observed in elementary collisions can be quantitatively explained 
[T7] as an effect of the different strange quark mass. 

These attempts relating the observed thermal features in hadron production pro- 



25 



cesses to quantum chaos or Hawking- Unruh radiation are still in a developmental stage. 
Whether they will keep their promises will be seen in the future. Certainly, both share 
the vision that there is a fundamental quantum mechanical mechanism behind this 
phenomenon and no (or little) room for a classical collisional thermalization process. 

11 Summary 

The statistical model has been applied to a wide variety of small and large systems, 
spanning more than two orders of magnitude in center of mass energy yfs and reveal- 
ing intriguing universal features of the hadronization process. The assumptions of the 
model were found to hold in a remarkable way for relative abundances and transverse 
momentum spectra of both light and heavy flavored species. 

This model is based on few simple principles which seem to capture key universal 
features in the process of hadron formation. While this does not answer the central 
questions of confinement and chiral symmetry breaking, which are of fundamental im- 
portance for QCD, it certainly allows us to understand some of the key properties of 
QCD in the non-perturbative regime. 
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